Numerical Renormalization Group at marginal spectral density: 
application to quantum tunneling in Luttinger liquids 
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Many quantum mechanical problems (such as dissipative phase fluctuations in metallic and 
superconducting nanocircuits, or impurity scattering in Luttinger liquids) involve a continuum of 
bosonic modes with a marginal spectral density diverging as the inverse of energy. We construct 
a Numerical Renormalization Group in this singular case, with a manageable violation of scale 
separation at high energy, capturing reliably the low energy physics. The method is demonstrated 
by a non-perturbative solution over several energy decades for the dynamical conductance of a 
Luttinger liquid with a single static defect. 



Bosonic description of fermionic systems, possibly sub- 
ject to strong interactions, has a long history, ranging 
from phase fluctuations in superconducting circuitsii^ fol- 
lowing Josephson's initial ideas, to quantum transport 
in metallic grains^ and in strongly correlated materials 
near the Mott transition^ where the phase conjugate 
to the electron charge is the relevant physical variable 
to understand the interplay of tunneling and Coulomb 
blockade. Another interesting example concerns one- 
dimensional electronic wires, the so-called Luttinger liq- 
uids (LL), where non- interacting plasmon modes pro- 
vide a faithful representation of electronic density fluc- 
tuations]^"— Quite remarkably, all these different physi- 
cal problems share very common features, because one 
can describe the electrons by a bosonic variable $ con- 
jugate to the electronic charge transfered in a nanoscale 
junction (or also to the charge density in a wire), via a 
phase factor e**. This mapping can be used for instance 
to represent the Josephson current in superconducting 
junctions (or also the fermionic fields in the bosoniza- 
tion language). This implies in turn that the nature of 
the phase dynamics determines the underlying physics: 
wild fluctuations of $ occur for instance in the presence 
of strong Coulomb blockade, leading to a rapid decay 
of the phase, and implying electronic localizatioui^i^ In 
contrast, phase localization is associated to small fluctua- 
tions of the electron charge, characterizing dissipationless 
supercurrenti or Fermi liquid statesi^ The intermediate 
situation of soft (algebraic) phase decay leads to the well- 
known non-Fermi liquid features of a LL»^"— 

In all these situations, great complexity arises due to 
the coupling of the bosonic mode to static disorder or 



ence of exponentially small energy scales arising from the 
impurity fiii^ for which powerful numerical methods have 
been developed in the past. The Density Matrix Renor- 
malization Groupie can tackle correlated wires, but on a 
linear energy scale only, which does not allow to really 
extract critical exponentsj^ the Numerical Renormaliza- 
tion Group (NRG) ^^d^ can however deal with impurity 
physics on an exponential energy range, but only for un- 
correlated Fermi liquids. A method that could incorpo- 
rate both virtues would therefore be quite useful, which 
is the goal of this Letter. 

Using the bosonic language, the description of inter- 
acting electrons by non-interacting bosons helps tremen- 
dously, but difficulties still arise. Apart from perturba- 
tive analysis (or fine tuning of the model parameters to 
allow an exact solution)}^"— the analytical bosonization 
technique offers limited information on quantum impu- 
rity problems, because the physics crosses over from weak 
to strong coupling, for instance due to Kondo screen- 
ing. Actually, all this complexity is already encoded by 
a static defect in LL, which drives the conductance from 

/h to zero on an energy scale that can be exponentially 
small in the backscattering amplitude, a problem that 
has triggered substantial work, based on approximate an- 
alytical method a^^d^'^"" — or numerical techniques on a 
linear energy scale, such as quantum Monte Carlo<^i^ 

The idea we henceforth present here is to use recent 
developments of the bosonic NRG^"— in order to tackle 
numerically the phase fluctuation problem in a broad 
range of parameters, with possible extensions to dy- 
namical defects. This however faces an immediate and 
seemingly intractable difficulty. Quite common to the 



dynamical defects, such as discrete Andreev levels^ in Josephson effect in a dissipative environment to quan- 



superconducting weak links^^) or magnetic Kondo im- 
purities in metallic junction o^^d^ and interacting unidi- 
mensional wireSi ^^d^ Focusing the discussion on the case 
of impurity effects in LL, but keeping this more general 
framework in mind, many technical and physical ques- 
tions are still open to date, both in the original fermionic 
formulation and in the bosonic version of the problem. 
On the fermionic side, one needs to handle strong interac- 
tions within unidimensional wires together with the pres- 



tum tunneling in resistive circuits) ^^i^^" — or to tunnel- 
ing into LLfi^ii^iiii^Sr— is the marginal form of the bare 
local bosonic spectrum, given by the correlation func- 
tion Q^^{iuj) = at imaginary frequency. The key step 
in the NRG procedure is the scale separation that re- 
sults from a logarithmic discretization of the energy band 
LUn = WcA~", with 1 < A. A generalized power-law den- 
sity of states (with exponent s > — 1 and high energy cut- 
off ujc) of the form J(a;) ~ 27ra;~^"'*ci;''8(a;c — uj) can be 
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considered both for fcrmionio^ and bosonic models^ — 
providing the fohowing eouphng strength of the states at 
energy w„: 
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27:^—^ A-"(^+i). (1) 
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For ah s > —1, the couphngs 7^ decay exponentiahy 
with n, which ahows an iterative diagonahzation of the 
problem: the possibihty of building progressively the 
Hilbcrt space from high to low energies is the reason be- 
hind the huge success of NRG to solve quantum impu- 
rity problems in a linear numerical effortJ^i^^ However, 
the marginal case s = — 1 is special in the sense that 
the couplings 7^ = 27rlog(A) do not decay anymore, in- 
validating clearly the whole scheme. We stress that we 
are considering quantum impurity Hamiltonians that de- 
pend explicitly on the phase factor e**, and not on the 
spatial derivative of the bosonic mode, dx^- This lat- 
ter case, which arises for instance in the so-called ohmic 
spin-boson model;^ corresponds to the much simpler sit- 
uation of linear spectrum (s = 1) of the field dx^, and 
can be easily handled by the bosonic NRG. 

This complete violation of scale separation for the 
marginal case s = — 1 seems to disqualify our proposed 
extension of the NRG. However, a free electron wire with 
constant density of states (described by the standard 
fermionic NRG at s = 0) is equivalent to a free bosonic 
bath with s = — 1 due to the bosonization mapping, so 
that one may believe that the marginal situation could be 
tackled using some clever variant of the bosonic NRG. In 
order to move forward, let us investigate with greater de- 
tail the problem of tunneling in LL. The fermionic Hamil- 
tonian reads in terms of second quantized left and right 
moving electron modes tp]^ fl(^) linear position x in 
the wire (omitting the electron spin for simplicity): 

H = J dxlivpipidxipL - ivFipRdxipj^ + g2tpRtpR'4'l'>pL\ 

-y6,[?/']jV'L+4^i?J|^=0 (2) 

where vp is the Fermi velocity, (72 the short-range 
Coulomb repulsion between left and right moving elec- 
trons, and Vbs the impurity backward scattering ampli- 
tude at the X = location of the defect (forward im- 
purity scattering and 174 interaction within a given Fermi 
point do not affect the physics, and were discarded.). The 
presence of the interaction term g2 clearly prevents a di- 
rect fermionic NRG solution of the model, which requires 
Fermi liquid leads. Yet, one can use the exact bosoniza- 
tion mapping^"— to re-express the electronic variables in 
terms of non-interacting collective charge density exci- 
tations $(a;) and conjugate field H(x). After standard 
manipulations^^— one obtains 

/fir 
— {[H(a;)]2 + [9,$(x)]2} ~v:cos[Vk^{x = 0)]t 

(3) 



where normal ordering of the cosine operator, which will 
be crucial for the rigorous formulation of the NRG algo- 
rithm, has been emphasized. We have also introduced 
a small backscattering energy scale v oc Vbs and the im- 
portant Luttinger liquid parameter K = [(1 — (72)/(l + 
.92)]^^^ < 1, into which all interaction effects have been 
encapsulated. 

Let us now present how the bosonic NRG^ can be 
tailored to address the impurity model (jSj, which has 
the form of a boundary Sine Gordon Hamiltonian. The 
derivation of the "star" -NRG follows the usual proce- 
dur e^^'^^ by considering the equivalent energy represen- 
tation in terms of a continuum of canonical bosons a| : 

H= de e ala^ - vtcos[y/K^]t, (4) 
Jq 

$ = $(a; = 0) = \/2 Tde 4±^. (5) 
^0 Ve 



The bosonic fields are then decomposed in Fourier modes 
(p € Z, n S N) on each interval ujn+i < e < a;„ of width 
dn = (l-A-i)A-": 



(6) 



The first NRG approximation consists in neglecting all 
p 7^ modes, keeping only the operators aj^ = aj^ q {^^^^ 
step becomes exact in the A — >■ 1 limit^i^). This leads to 
the "star" -Hamiltonian: 
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: (7) 



with the "impurity" coupling strength already given in 
Eq. (IlJ by 7,^^ = 27rlog(A) in the marginal case s = — 1. 
The typical energy ^„ in each shell is defined by: 



1 /■"" 1 - K-^ 



.A- 



(8) 



As a benchmark of the discretization for the marginal 
case s = — 1, one can easily compute from ([7|) the result- 
ing approximation for the original Green's function: 



+00 



WrA 



n=0 CJ-^ + 
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(9) 



which can be checked to converge exponentially fast at 
cj ^ to the exact result Q%{ioj) = even for A = 2 
(we keep this standard value from now on). However, de- 
spite the clear exponential decay of the energies (|S|) , the 
non-decreasing value of the couplings 7„ implies a vio- 
lation of scale separation on all shells, and prevents the 
solution by iterative diagonahzation of Hamiltonian ([7]) . 

The first key idea in successfully constructing the 
marginal bosonic NRG is to assume that the energy spec- 
trum is also bounded from below: 
2tt 

J(w) = — 6(i:Jc - w)B(a; - Winin)- (10) 
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Clearly both the energies e„ and the couplings 7„ are not 
modified by this choice (7„ still do not decay), and they 
arc just cut off for n > riniin, with ~ WcA~"™'", 

so that nothing seems gained naively. We can however 
try to pursue with the second step of the standard NRG 
procedure, which amounts to the exact mapping on the 
Wilson chainiiS— This simple tridiagonalization proce- 
dure of Hamiltonian ([7]) leads to the following "chain" 
form in terms of new canonical bosons 61 : 



+ 00 

n=0 



-vlcos 



[bl + bo) 



(11) 



with the parameter tjq = 27rlog(wc/'^min)- Clearly, the 
impurity part of the chain Hamiltonian (jlip breaks down 
for Wmin 0, Owing to the divergence of ijo, but one can 
check numerically that the construction is valid for non- 
zero ojmin. The on-site energies e„ and hoppings t„ of the 
Wilson chain can indeed be obtained by numerical tridi- 
agonalization of Eq. (O. For the value Wmin = 10~^ of the 
lower cutoff, these are plotted together with the star pa- 
rameters on Fig. [T] The exponential decay of both chain 
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FIG. 1: (color online) Left panel: parameters ^„ and 7„ of 
the star-NRG as a function of n for < n < nmin ~ 16; the 
coupling 7„ does not decay and violates scale separation on all 
shells. Right panel: parameters e„ and t„ of the chain-NRG; 
scale separation is only broken on the first shell, as seen by the 
initial increase of both parameters, before further exponential 
decay (shown by dotted lines as guides to the eye). 

parameters e„ and t„ that we discover here is clearly a re- 
markable surprise, that enables the extension of the NRG 
to the marginal situation s = — 1. This crucial feature 
comes at a small price, seen by the first increase of the 
chain parameters from site n = to site n = 1. Thus the 
maximal violation of scale separation in the star NRG 
presents a small remanence in the chain NRG, limited 
only to the first shell. Interestingly, the initial jump of 
the parameters is just proportional to log(ti;c/^i^min), so 
that the lower cutoff Wmin can be decreased on exponen- 
tial scales without paying a huge numerical cost. 

A last difficulty due to the unusual form of the im- 
purity Hamiltonian (jlip must be addressed. In the 
standard NRG ;^^'^^ only linear to quadrilinear opera- 
tors are present in the Hamiltonian. However, the cen- 
tral physical role played by the phase factor e'* leads 
to a cosine term at the impurity site, hence to an oper- 
ator of infinite order, which by the bosonization rules^ 



must also be normal ordered. For a generic operator 
O — ; cos[a(5j-l-6o)] * this reads O = cos(q;6J) cos(a6g) — 
sin(a6j) sin(a6Q). Using the Fock states |m) of the 
bosonic creation operator 6j on the initial site n = 
of the Wilson chain, one obtains the matrix elements: 







Min(m,p) 

(m|0|p) = i/mlp! Tie ^ 



(ia) 



m-\-p—2k 



k=0 



(m - k)\{p - ky.kl 



(12) 



The construction of the impurity term in proceeds 
by a truncation of the infinite Fock space on the initial 
Wilson site limited to states with occupation number less 
than a given iVoi s-nd use of the matrix elements (jl2p . 
Typically iVg = 150 ensures a good representation of the 
Hamiltonian. Each further site ti > of the chain is de- 
scribed by a basis of Nj, states^^ (we take Ni, = 12 here). 
At increasing n, the growing size of the total Hilbert 
space becomes rapidly unmanageable, and a truncation 
to A^trunc states is required (this approximation is com- 
mon to all NRG schemesiS). Typically A^trunc = 800 was 
employed in all further computations, and we also set 
= 1 as the basic energy unit. 




FIG. 2: (color online). Bosonic correlation function Re[t/#(cj)] 
at real frequency u for the LL parameter K = 1/2 comparing 
(bottom to top) the NRG to the exact result (|13p . and to 
the strong and weak interaction perturbation theory given 
respectively by (|14p and (|15|) (these two expressions are by 
accident equivalent for K — 1/2, but nonetheless not exact). 

In contrast to more complex extensions of impurity 
models with dynamical degrees of freedom (such as the 
Kondo model in a Luttinger liqui d^'^d^ ), the present im- 
purity problem benefits from several known limits, that 
allow to benchmark our numerics. For instance there ex- 
ists an exact solution for the dynamical conductance^"— 
(in units of e'^/h) at K = 1/2: 

'(c.) = ^UeiG^iu:)] = i - ^atan (^) (13) 

where fi = e'^v'^, with Euler's constant 7, is the crossover 
energy at which the impurity cuts the chain (for K — 
1/2). Perturbation theory works also at strong interac- 
tion K <ti 1, in which case the self-consistent harmonic 
approximation appliesi^ 
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FIG. 3: (color online). Dynamical conductance G(a;) 
in units of e^/h for several values of the LL parameter 
K = 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9 obtained by the 
NRG (symbols, bottom to top). Comparison is made for 
K < 0.5 to the strong interaction and for K > 0.5 to the weak 
interaction limits. Perturbative methods cannot describe ac- 
curately the NRG solution beyond their applicability range. 



with O* = 271 Kv i-A' the crossover scale. Finally, the 
limit of weak interaction 1 — if ^ 1 is also known from 
several approaches j^ii^iiiSi"— 



KujK 



(15) 



^*-2 + (r!-)*-2 

Fig. [2] compares om- NRG data for if = 1/2 with the ex- 



act solution and perturbation theory, which shows the ex- 
cellent convergence of the NRG and the sizeable discrep- 
ancies of both perturbative expansions. More systematic 
analysis for various K values in Fig. [3] demonstrates the 
progressive departure of the perturbative results from the 
numerical solution. The ability of the marginal bosonic 
NRG to describe non-perturbatively universal transport 
features with high accuracy should therefore make it a 
precious tool to test scaling behavior of impurity physics 
in LL. 

To conclude, we have established an extension of the 
NRG to deal with the marginal situation of a density 
of states diverging as the inverse of energy. The po- 
tentially most promising applications of the NRG at 
marginal coupling concern the physics of dynamical im- 
purities coupled to phase fluctuations, a large class of 
physical problems where no alternative analytical or nu- 
merical techniques exist to date. This development could 
allow to address many currently open issues, such as non- 
equilibrium transport with strong correlations (using a 
mapping onto equilibrium q-oscillator models^), Kondo 
physics in Luttinger liquidsf^^— and ohmic dissipation 
in Andreev level qubits4"— 
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